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Abstract. Closed-orbit theory provides a general approach to the semiclassical 
description of photo-absorption spectra of arbitrary atoms in external fields, 
the simplest of which is the hydrogen atom in an electric field. Yet, despite 
its apparent simplicity, a semiclassical quantization of this system by means of 
closed-orbit theory has not been achieved so far. It is the aim of this paper 
to close that gap. We first present a detailed analytic study of the closed 
classical orbits and their bifurcations. We then derive a simple form of the 
uniform semiclassical approximation for the bifurcations that is suitable for an 
inclusion into a closed-orbit summation. By means of a generalized version of 
the semiclassical quantization by harmonic inversion, we succeed in calculating 
high-quality semiclassical spectra for the hydrogen atom in an electric field. 
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1. Introduction 

Because it is classically integrable, the hydrogen atom in a homogeneous electric field 
can be regarded as the simplest and most fundamental atomic system apart from 
the field-free hydrogen atom, and consequently the semiclassical treatment of that 
system is of particular importance. Due to the integrability, of course, the Einstein- 
Brillouin-Keller torus quantization rules apply, and a semiclassical calculation of 
photo-absorption spectra based on the existence of quantized classical tori was indeed 
carried out by Kondratovich and Delos 

A more general approach to the semiclassical description of atomic spectra is 
furnished by closed-orbit theory |^,|^], which is intended to be applicable to atomic 
systems exhibiting either regular, chaotic or mixed classical behaviour. Therefore, it 
is of vital interest to see how it can be applied to this apparently simple example. 

Closed-orbit theory provides a semiclassical approximation to the quantum 
response function 



PACS numbers: 32.60.+i,03.65.Sq,31.15.Gy,32.70.Cs 
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which involves the energy eigenvalues i?„ and the dipole matrix elements {i\D\n) 
connecting the excited states n to the inital state i. Semiclassically, the response 
function splits into a smooth part and an oscillatory part of the form 

5°^^(i?) = ^A.o.(i?)e'^-(^) , (2) 



CO. 



where the sum extends over all classical closed orbits starting from the nucleus and 
returning to it after having been deflected by the external fields. Sc.o. is the classical 
action of the closed orbit. The precise form of the recurrence amplitude .Ac.o. depends 
on the geometry of the external field configuration |§, ||, ^ . In the presence of a 
single external field, it reads 

Ao=4.^^:M4MeUV.). (3) 
\mi2\ 

if the orbit is directed along the field axis and 



A.O. =2(2.)3/^^!^^5r(^,)y(^.)exp(if. + iJ) (4) 

if it is not. Here, and d / are the inital and final angles of the trajectory with respect 
to the field axis, fi is its Maslov index, and mi2 an element of the monodromy matrix. 
The function 3^(??) encodes the relevant properties of the inital state and the exciting 
photon. The amplitudes are independent of whether the classical dynamics is regular, 
chaotic, or mixed. This is in stark contrast to semiclassical trace formulae 0|, which 
represent the quantum density of states as a sum over classical periodic orbits formally 
analogous to (H), but require different semiclassical amplitudes in the cases of regular 
or chaotic dynamics. This property makes closed-orbit theory an ideal starting point 
for a semiclassical quantization of typical atomic systems in external fields. 

Nevertheless, even for the simplest of such systems, the hydrogen atom in 
an electric field, any attempt to semiclassically quantize based on closed-orbit 
theory so far has been unsuccessful. The reason for this failure is that, even 
though the closed orbits of the Stark system are easy to describe and classify 
[H, a multitude of bifurcations exist. Bifurcations of closed orbits lead to a 
divergence of the semiclassical recurrence amplitudes and thus spoil the semiclassical 
spectrum. Uniform semiclassical approximations 11| smooth these divergences. 



but they have so far been restricted to the calculation of low-resolution spectra 
and have prevented the use of any high-resolution semiclassical quantization scheme. 
Only recently were we able to introduce a generalization of the semiclassical 
quantization by harmonic inversion which is capable of coping with uniform 
approximations. We have chosen the hydrogen atom in an electric field to illustrate 
our quantization scheme because its closed orbits and their bifurcations can easily be 
understood and its simplicity makes it a archetypal problem of closed-orbit theory. On 
the other hand, bifurcations of closed or periodic orbits are a generic phenomenon of 
mixed regular-chaotic systems, so that their treatment in a semiclassical quantization 
scheme is an important task with widc-rangig applications. The principal difficulties 
of this task are already present in the Stark system, so that it furnishes an ideal testing 
ground for the novel technique. 

In we focused on the presentation of the quantization scheme itself. The 
central theme of the present paper is the classical and semiclassical dynamics of the 
hydrogen atom in an electric field. In section ^ the classical dynamics of the Stark 
system is described and the equations of motion are solved explicitly. In section ||, 
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the closed orbits and their bifurcations are discussed. Semi-analytic formulae for all 
relevant orbital parameters are derived which have not been given in the literature 
so far. In section ^ a uniform approximation for the closed-orbit bifurcations in the 
Stark system is constructed. Although this has been done before we derive 

our result in a much simpler way as was done previously and bring it into a form 
that is considerably easier to apply. Section |^ concludes by presenting both low- and 
high-resolution semiclassical spectra. They prove that uniform approximations can 
indeed be included into a high-resolution semiclassical quantization. These results 
solve the long-standing problem of a closed-orbit theory quantization of the Stark 
effect and serve as a prototype example for the application of the novel method to 
other systems. 



2. The equations of motion 



If we assume the electric field F to be directed along the z-axis, the Hamiltonian 
describing the motion of the atomic electron reads, in atomic units, 

H=Ip^^- + Fz , (5) 
2 r 

where = -\- -V . It can be simplified by means of its scaling properties: 
If the coordinates and momenta are scaled according to i, — F^^'^x^ p — F~^/'^p, 
the classical dynamics is found not to depend on the energy E and the electric field 
strength F separately, but only on the scaled energy E — F~^/'^E. In this work, 
we will use scaled classical quantities throughout in sections |^ and where purely 
classical calculations are presented. For the sake of simplicity, we will not mark them 
with a tilde. From section |^ on, the notational distinction between scaled and unsealed 
quantities will be taken up again. 

As is well known the Hamilton- Jacobi equation corresponding to (||) 

separates in semiparabolic coordinates pi,p2, where 



pi = \Jr ~ z , Pi = \Jr + z , (6) 

and Lp is the azimuth angle around the field axis. Notice that p\ and pi are defined 
by (H) up to a choice of sign only. 

The Hamilton- Jacobi approach leads most naturally to the computation of the 
action-angle variables fundamental to the EBK quantization. For a discussion of closed 
orbits, however, it is essential to follow the temporal evolution of a trajectory. This can 
most conveniently be achieved if the trajectories are parameterized by a parameter t 
related to the time t by 

dT^2rdt. (7) 

In the following, primes will be used to denote differentiation with respect to r. The 
Hamiltonian describing the dynamics as a function of r is given by 

n^2r (H-E) 

= \ {pI. + pI + +Pl)-pt + pt-i)=0. 

(See jl^ or j|] for a discussion of how to treat the pseudotime transformation (0) 
within the framework of Hamiltonian dynamics.) 
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Separating the pseudotime-Hamiltonian (|^) into pi-dependent and p2-dependent 
parts and using Hamilton's equation of motion p'^ — dTi./dpp- — Pp., we obtain 

2 

= ±p| + 2i?p2 _ ^ + c, for J = 1,2 (9) 

with separation constants ci and C2 related by 

ci + C2 = 4 . (10) 
Equation can be integrated to yield 

dpi _ f dp2 ^^^^ 



pi + 2Epl-^+ci ^-pi + 2Epl-^ + C2 

which gives the parameter t as an elliptic integral in the coordinates pi and p2. 

In the case of closed orbits, = 0. At the nucleus, we have p'^ = Cj by (||). The 
initial conditions for an orbit starting at the nucleus at time r = thus read 

pi(0) = , p2(0) = , 

p'i(o) = v5r, p'2{o)^v^. 

The constants ci and C2 describe the distribution of energy between the uphill and 
the downhill motion. 

According to (^), for real orbits pj can only assume values which make 

f,ip,)^±p^ + 2Ep'j+c,>Q. (13) 

For the motion in the pj-direction to be bounded, there must be a real turning 
point, where fj{pj) — 0, because otherwise pj will keep increasing forever. From 
the factorized form 

with 



(12) 



pI± ^~E± ^/E^ - ci , (15) 

it is obvious that for _E < the zeros of /i(pi) are real if 

ci < E^ (16) 

and complex otherwise. Therefore, a closed orbit can only exist if (|l^ ) holds. For 
energies E below the Stark saddle energy Es — ~2, this condition is satisfied for all 
real orbits because C2 = 4 — ci > 0. For E > Ec,, there must be a sufficiently strong 
component of the motion in the uphill direction so that the energy of the downhill 
motion does not suffice to cross the saddle. 
For the p2 motion we obtain similarly 

f2{p2)^-{pl-pl+){pl-pl-) , (17) 

with 

pI^^E± + C2 . (18) 

The coordinate p2+ is real for any C2 > 0. Therefore, the p2-niotion is always bounded 
for real orbits. 
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With the help of the factorizations ( |14| ) and (|l^), the elliptic integrals (|ll|) can 
easily be reduced to Legendre's standard integral of the first kind 

1 



p+ V 

The parameters are 



P+ 

arcsm I pi 

P+ V VVci 



pI 



mi = 

and the abbreviations 

P+ = Pi 



mi 



ci 



P- 



■ T ( arcsin ( ^ p2 



m2 



9 4 



m2 = ^— 
P2- 



C2 

„4 



P- 



-Pt 



(19) 



(20) 



(21) 



were introduced. Equation (|19|) can be solved for pi and p2 in terms of Jacobi's elliptic 
functions [|6| to yield 



(22) 



Pi('^) = -^sn (p+t| mi) , 
P+ 



P2{t) 



P- 



sn (p-t| m2) 



These results incorporate the initial conditions (O 



3. Closed orbits 

Closed orbits are characterized by the condition that at a pseudotime tq > the 
electron returns to the nucleus, so that 

Pi(to) = P2(to) = 0. (23) 

In the simplest cases, either pi or p2 vanishes identically. The electron then moves 
along the electric field axis. 

If ci = 0, the downhill coordinate pi is zero. The electron leaves the nucleus in 
the uphill direction, i.e. in the direction of the electric field, until it is turned around 
by the joint action of the Coulomb and external fields and returns to the nucleus at a 
pseudotime 

/C(m2) 



T2 



(24) 



where /C(m) denotes the complete elliptic integral of the first kind |T^. Thus, T2 
corresponds to a half period of p2 and, due to the freedom in the choice of sign, 
to a full period of the uphill orbit in position space. The uphill orbit is repeated 
periodically and closes again at tq = Iti with integer I. 

The second case of axial motion is obtained if ci =4, which corresponds to a 
downhill motion opposite to the direction of the electric field. The pi motion closes 
at pseudotimes t = fcri with 

/C(mi) 



P+ 



(25) 



if the energy E is less than the saddle point energy E'^ — —2. li E > Es, the electron 
crosses the Stark saddle and the orbit does not close. 
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In the case of a non-axial orbit, the uphiU and downhill motions must close at 
the same time 

To = fcri = It2 . (26) 

If ( p6| ) is satisfied for any given integer values of k and I, the orbit returns to 
the nucleus after k half periods of pi and I half periods of p2, corresponding, in 
Cartesian coordinates, to k full periods in the downhill and I full periods in the uphill 
direction. We therefore refer to k and / as the downhill and uphill repetition numbers, 
respectively, and identify a non-axial closed orbit by its repetition numbers (fc, I). 

If the scaled energy E is fixed, the half-periods ri and T2 depend on the separation 
constant ci. For given k and (^6|) can therefore be read as an equation for ci, thus 
determining the initial conditions of a closed orbit. We are now going to investigate 
the solution of this equation to determine the conditions for a (fc, I) orbit to exist. 

Figure |l| shows the dependence of ri and T2 on ci for a fixed energy E. For all ci , 
the uphill period T2 is smaller than the downhill period ti. Furthermore, as ci — > —00 
we obtain from (§3) and (|2^), to leading order. 



whereas 



^ 2/C(-l) 

-i--^-(Z^' (27) 
^ 2/C(-l) 

-1 ~ (Z^ ' (28) 

so that the ratio ti/t2 approaches 1 from above. On the other hand, ri diverges at 
Ci = E'^, so that ti/t2 ^ 00 as ci ^ E^ . Thus, for any energy E < and positive 
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Table 1. Comparison of parameter ranges for real and ghost orbits 



integers fc, I with I > k, there is a unique solution ci to (|2^). It is in the range 
— oo < ci < E^. There is no solution for / < k, since ri is complex for ci > 
whereas T2 remains real. Therefore, orbits (fc, I) with Z < fc do not exist. 
A real closed orbit must satisfy 

< Cl < 4 (29) 



due to (^^ and (y_2|). If (29) is not met, the initial conditions are complex and a 
ghost orbit in the complexified phase space is obtained. More precisely, if ci > 4, the 
initial velocity /92(0) is imaginary, hence P2{t) will be imaginary at all times. Similarly, 
Pi{t) will be imaginary if ci < 0. In both cases, if the complex conjugate of the ghost 
orbit is taken, one of the semiparabolic coordinates changes sign. This change of 
sign does not alter the Cartesian coordinates, and consequently all ghost orbits are 
invariant with respect to complex conjugation. By the same token, the starting angle 
of a ghost orbit is defined up to complex conjugation only. We will always choose 
the imaginary part of starting angles to be positive. The properties of real and ghost 
orbits are summarized in table |l]. 

To investigate the bifurcations of the closed orbits, we will now discuss the 
dependence of ci on E for fixed repetition numbers fc, I. For any energy E and ci = 0, 
the downhill period ri is equal to the period 

rcoui = ^ (30) 

of the pure Coulomb dynamics. For ci = 4, the uphill period T2 equals tcouI- For 
any other value of ci, both ti and T2 converge to rcoui as E — oo, whence in this 
limit ri/r2 — > 1. Thus, ti/t2 < ?/fc at ci = 4 if i? is sufficiently low. The solution 
to ( p6| ) must then lie in the interval ci € [4, i?^], so that the (fc, ?)-orbit is a ghost. 
As E increases, ci decreases. At the critical energy Eg^n where ci = 4, a real (fc, I) 
orbit is generated. It bifurcates off the downhill orbit that is located at ci = 4. For 
E > -Egen, Cl decreases further. When the energy i?dest where ci = is reached, the 
(fc, I) orbit collides with the uphill orbit and becomes a ghost again. As the singularity 
of Ti approaches zero a.s E y 0, there is an Edest < for any (fc, I). 

The bifurcation energies i?gen and £^dest can be determined from ( p6| ) if ci = 4 
is prescribed for the generation of an orbit, or ci = for its destruction. For these 
cases, equation (|26|) simplifies to 



I 2'V2 
fc ' 




(31) 



in terms of the dimensionless variable e = 4/£'^ = (Es/E)'^ for the generation of the 
{k,l) orbit and 
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Figure 2. Orbital parameters close to the bifurcations of the orbit (4,5). 
Solid line: non-axial orbit (4,5), long-dashed line: fifth uphill orbit, short- 
dashed line: fourth downhill orbit. Vertical lines indicate the bifurcation energies 
E„cn = -2.3597 and 

^dcst — — 1.3790, respectively 



for its destruction. These equations provide a simple and stable method to determine 
the bifurcation energies. They can be expected to yield more accurate results then 
the numerically computed monodromy matrix elements used by Gao and Delos 

The sequence of bifurcations described above is illustrated in figure where 
characteristic data of the orbit (4, 5), and the downhill and uphill orbits it bifurcates 
from, is shown. The bifurcation energies -Egcn = —2.3597 and i?dcst = —1.3790 are 
characterized by the conditions Ci = 4 and ci — 0, respectively. The transition 
between real and ghost orbits can most clearly be seen from the energy dependence of 
the starting angle "di , which is real between the bifurcation energies and acquires a non- 
zero imaginary part outside this interval. In addition, at the bifurcation energies the 
actions and orbital periods of the non-axial orbit coincide with those of the appropriate 
axial orbits. Note that the action of the non-axial orbit is always smaller than that of 
the axial orbits. 

It is also apparent from figure ^that the monodromy matrix element mi2 vanishes 
at the bifurcations. For the axial orbits, further zeros of mi2 appear in figure ^. They 
correspond to additional bifurcations these orbits undergo. As fc < ^ for a non-axial 
orbit, the Ith repetition of the uphill orbit undergoes / — I bifurcations, where the orbits 
(fc, I) with k = 1, ... ,1 — 1 are destroyed. All four bifurcation energies of the fifth uphill 
orbit are visible in figure |[ A downhill orbit, on the contrary, undergoes infinitely 
many bifurcations before it ceases to exist at E = —2. Only the first members of this 
infinite cascade are resolved in figure 0. 
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The separation parameter ci for a closed orbit must be calculated numerically 
from (p6|). Once it is known, all orbital parameters can be found analytically. The 
pertinent formulae will now be derived. 

First of all, if we transform the orbits ( p^ ) to Cartesian coordinates, we find that 
the starting and returning angle ■& of a closed orbit, measured with respect to the field 
axis, is related to the separation constants by 

2cos- = Vci, 2sm- = Vc2, (33) 

so that 



d = 2 arccos ^ . (34) 

The angle determined from (|3j) is obviously real and confined to the interval 
< 1? < TT if < ci < 4, as was to be expected for a real orbit. For a pre-bifurcation 
ghost orbit, i.e. ci > 4, the equations ( |3^ ) can be satisfied if d = ia is chosen purely 
imaginary with 

2 cosh — = ^/cT , 2i sinh — = ^/c^ . (35) 

This choice of "d, which is unique up to the addition of multiples of Air, makes 
continuous at the bifurcation. Similarly, for ci < 0, i.e. for a post-bifurcation ghost 
orbit, we set = n + ia with 

2isinh— = 2cosh— = -^€2. (36) 

Again, the real part makes ?9 continuous at the bifurcation. In both cases we assume 
a > 0, which is possible because the orbit is invariant under complex conjugation. 
Using pp^ — p'j , we obtain the action integral 



S ^ J Ppi dpi + J pp^ dp2 

= 2k p[ dpi +21 p'2 dp2 

Jo Jo 



(37) 



The equations of motion (0) then yield 



\/cT/p+ rV^/p 



S^2k ^fi{pi)dpi+2l / VhiP2)dp2 

Jo Jo (38) 

2 2 
^ -kpll{mi)- -lpmm2) 

with 

/(m) = (m- 1)/C(m) + (m + l)5(m) (39) 

and the complete elliptic integrals of the first and second kinds JC{m) and £(m) |p^ . 
Note that /(m) < if m < 0, so that both terms in (|3^ ) are positive for real orbits. 
For a downhill orbit, the uphill repetition number I is undefined. As, however, in this 
case m2 = and /(O) = 0, the uphill contribution drops out of ( |3^ ) for a downhill 
orbit. Similarly, the downhill contribution vanishes for uphill orbits. The action S 
given by (^8|) is real for both real and ghost orbits. 

The physical-time period of a closed orbit is given by 



r« 2 r" 2 

T = Pi dpi + P2 dp2 
Jo Jo 



(40) 



2kp+ (/C(mi) — £{mi)) — 2lp^ (/C(m2) — £(7712)) 
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As is the action, the period T is always reaL 

The calculation of the monodromy matrix element is more difficult . 
proceeds from the 2x2 Jacobian matrix 

dp{To) 



J 



dp{0) 



It 



(41) 



that describes the change of the final position upon a variation of the initial 
momentum. The initial and final momenta are given by 



p(0) 



p(to) 



(-1)' 



(42) 



They have norm 2 by (|1C|). As the monodromy matrix characterizes variations 
perpendicular to the orbit, toi2 is obtained from J by projecting onto unit vectors 
perpendicular to the momenta (^2|). Thus, 



mi2 = i ( (- 



-1)'V^, -(-1)^=751 )• J- i 



C2 
Cl 



Since the dynamics is separable, the matrix J is diagonal, whence 



, -.sfcCi dp2{T) 

^ ^4 dp2{0) 



For the remaining derivatives we find from (p3) 



dpi{0) \dpiiO) p 



d 



Cl 



sn (p+t| mi) H 



d 



■ sn { p+t\ mi) 



(43) 



(44) 



(45) 



p+ dpi{0) 

A similar expression holds for p2- If a non-axial orbit {k,l) is considered, the first 
term in ( ^5|) vanishes at r = tq by virtue of the resonance condition (26). The second 
term can be evaluated in an elementary, but lengthy calculation to yield 

kmiC2 



mi2 



(-1) 
+ (-1) 



k+l 



Ap+VW^ 

k+l ^"l2Ci 



Cl 



with 



d(m) 



Ap^^/W^ 
E{m) /C(m) 



C2 



{2Ed{mi) - p\K{mi)) 
{2Ed{m2) - pllC{m2)) 



(46) 



(47) 



m(l — m) 

Prom ( ^6|) it can be verified that, up to a choice of sign, the matrix element mi2 
for the orbit {nk,nl), which is the nth repetition of {k,l), equals n times the matrix 
element for the orbit (fc, I), as has been shown previously by Gao and Delos Q by an 
abstract argument using the neutral stability of the orbits. It is also clear from ( ^6| ) 
that mi2 vanishes when the orbit undergoes a bifurcation, because mi — > as ci — )■ 0, 
and /C(mi) and d{mi) both approach finite limits. 

For an uphill orbit the second term in (45) vanishes because ci = 0, whereas 
the first term is nonzero in general because the axial orbits do not obey a resonance 
condition akin to (26). The derivative reads 

dpi{T) 



9^1 (0) 



ci=0 



sin(V^2i^T) 



(48) 
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so that 



(49) 

Similarly, the matrix element for a downhill orbit is 



,,fc sin(V-2^To) 

"^'^ = ^-'^ ^ • ^^'^^ 
Finally, the Maslov indices of the closed orbits need to be determined. For the 
uphill orbits, we proceed as follows: First, there is a caustic whenever the orbit reaches 
either the nucleus or the turning point, totalling to 2Z — 1 caustics. Second, the Maslov 
index increases by two (corresponding to two independent directions transverse to the 
orbit) whenever the orbit is intersected by neighbouring trajectories, i.e. when ( ^8| ) 
vanishes. This occurs, by (^), after a pseudotime tcouI, so that the intersections 
contribute 2[To/rcoui] to the Maslov index, [x] denotes the integer part of x. Third, 
if the semiclassical amplitudes are written according to the conventions of |^,|^, the 
Maslov index must be increased by one, so that for an uphill orbit it finally reads 

''^^^^ (51) 



^JL = 2[l + 



_ TCoul 

Using the same reasoning, we find 

^^2(k+\^]) , (52) 

for a downhill orbit. As the downhill period ti is always larger than rcoui, the downhill 
Maslov index (|2|) is equal to 4A; for sufficiently low E. 

The Maslov index of a non-axial orbit involves the number of zeros of mi2 
encountered along the orbit, which is hard to find from (^). Instead, we exploit the 
observation that the change of Maslov indices in a bifurcation can be determined from 
the normal form describing the bifurcation. A normal form suitable for the bifurcations 
in the Stark effect will be given in section ^ It predicts that in a bifurcation the Maslov 
index of the axial orbit increases by two, which is consistent with (^ij) and (|5|), and 
the Maslov index of the non-axial orbit takes the intermediate value. A schematic 
drawing of the shortest orbits and their bifurcations in shown in figure ^ It implies 
that the Maslov index of a non-axial orbit is 

^ = 2(A:-|-0 - 1 ■ (53) 

Since along the orbit there are k—1 zeros of pi and l — l zeros of p2, each corresponding 
to an intersection with the electric field axis and thus contributing 1 to the Maslov 
index, there must be fc -I- ^ — 1 zeros of mi2 . This result can be confirmed numerically. 

For ghost orbits, Maslov indices cannot be determined by counting caustics. 
The correct phase of the ghost orbit contribution is fixed by demanding that the 
semiclassical amplitude be continuous across the bifurcations, as required by the 
uniform approximation to be developed below. 



4. The uniform approximation 

In the following sections we present the semiclassical treatment of the Stark effect. 
From now on, scaled quantities will again be marked with a tilde. 

It has become clear in section ^ that the apparently simple dynamics of the 
hydrogen atom in an electric field contains numerous bifurcations of closed orbits. 
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Figure 3. Periods of the shortest closed orbits, labelled by their Maslov indices, 
as a function of energy (schematic). Solid lines: non-axial orbits, short-dashed 
lines: downhill orbits, long-dashed lines: uphill orbits. 



They lead to zeros of the monodromy matrix elements myi of the bifurcating orbits 
and thus to the divergence of the pertinent closed-orbit amplitudes (^) , so that they 
spoil the semiclassical spectrum (see section ^) . To overcome this problem, a uniform 
semiclassical approximation giving the collective contribution of the bifurcating orbits 
must be derived. For the bifurcations occuring in the hydrogen atom in an electric 
field, this was done by Gao and Delos [||. These authors constructed a semiclassical 
wave function near a bifurcation and essentially redid the derivation of closed-orbit 
theory using that wave function. Their uniform approximation was extended by Shaw 
and Robicheaux [|l^,|ri[ to include ghost orbits. 

In the form given in 10, ll|, the uniform approximation cannot be calculated 
from the actions and recurrence amplitudes characterizing the bifurcating orbits, but 
requires the determination of additional parameters. It is therefore hard to apply if 
many closed-orbit bifurcations are to be dealt with. In this section we will present 
a novel approach to uniformizing the bifurcations in the Stark system that is not 
only simpler than previous derivations, but also yields the uniform approximation in 
a more convenient form. It will be shown that only the actions recurrence strengths 
of the isolated orbits are needed to compute a uniform approximation. This makes 
the uniform approximation as easy to use as the simple isolated-orbits formula. 

The crucial step in the construction of a uniform approximation is the choice of 
a suitable ansatz function whose stationary points describe the classical closed orbits. 
In the Stark system, the structure of the bifurcations is determined by the rotational 
symmetry of the system. Non-axial orbits occur in one-parameter families generated 
by the rotation of any particular member around the electric field axis, whereas the 
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axial orbit is isolated at all energies except at a bifurcation. The normal form must 
share these symmetry properties. 

We choose the normal form $ to be defined on a plane with Cartesian coordinates 
(x, y) or polar coordinates 

X = r cos 95 , y — r sin 1^9 . (54) 

Due to the rotational symmetry, <i> must be a function of the radial coordinate r 
only, independent of the angle ip. To be a smooth function of x and y, it must be 
a function of r^. Any such function has a stationary point at the origin x = y = Q 
that corresponds to the axial orbit because it is rotation invariant. Any stationary 
point at a non-zero radial coordinate occurs on a ring r = r,, and describes a family 
of non-axial orbits. A normal form suitable to describe the bifurcations occurring in 
the Stark system must therefore have stationary points on a ring with radius Tc that 
shrinks to zero as the bifurcation is approached. 

These requirements are fulfilled by the normal form 

^a{x,y) = \r^ -\ar'' (55) 

with a real parameter a. Apart from the trivial stationary point at the origin, ( p5|) 
has a ring of stationary points at r — y/a, which is real if a > and imaginary if 
a < 0. Thus, (|5^) describes a family of non- axial real orbits present for a > which 
then contracts onto an axial orbit and becomes a family of ghost orbits. 
For the uniform approximation we make the ansatz 

^{E) = I{a)&^--^^'> , (56) 

where S'ax(-E) denotes the action of the axial orbit and 

I (a) = f dxdyp(x2 -f y2)ei*o(^,y) ^ (57) 



with an arbitrary smooth rotationally symmetric function p(r^). The parameter a 
and the function p(r^) must now be determined such that, when the distance from 
the bifurcation gets large, the uniform approximation reproduces the isolated-orbits 
formula. The asymptotic behaviour of ( ^ ) can be found if ( ^7|) is evaluated in the 
stationary-phase approximation. The contribution of the stationary point at the origin 
is straight-forward to evaluate. It yields 

/(«)l.=o = T^^'W^"'^''" = ~ ' ^^^^ 

where 

f : a < 

vo = <^ 59 

is the number of negative eigenvalues of the Hessian matrix. From this number it is 
apparent that the Maslov index of the axial orbit increases by 2 in the bifurcation. 

The contribution of the non-axial orbits is present if a > 0. Due to the rotational 
symmetry, a straight-forward stationary-phase approximation to this contribution fails 
because the stationary points are not isolated. In polar coordinates, the integration 
over the angle ip yields the constant 2tt, and the stationary-phase approximation can 
be applied to the remaining integral over r. This yields 

/(a)|,i„ =2^^/^p(a)e--V^ (gO) 
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Note that the contribution of the non- axial orbit remains finite as a ^ 0. 
Thus, the uniform approximation (p6|) asymptotically equals 



vi;(ij) ^ ^ — p{0) e'^^^(^) + 2^V^p(a) e'(5=^(^)-«'/4) . (61) 
a 

This result should reproduce the isolated-orbits formula 

AUE) e'^^''*^) + AnoniE) e'^""(^) (62) 

in terms of the actions and recurrence amplitudes of the axial and non-axial orbits. 
This can be achieved if we demand 

■S'non = ■S'ax - , (63) 

Ax = - — P(0), Anon = 271 V^i p{a) . (64) 

a 

The normal form parameter 



±2v^S'ax - 5'noii (65) 



can be determined from (|63|). The sign of a has to be chosen according to whether 
the non-axial orbits are real or complex. 

In the limit a — + 0, i.e. close to the bifurcation, p{a) « p(0). Equation (|^ thus 
imposes the constraint 

: Aion = -aAx (66) 



V27r 

on the semiclassical amplitudes. Note that both sides of this equation arc finite as 
o ^ 0. In particular. Anon is continuous at a = 0. This condition was used in 
section || to fix the Maslov indices for the ghost orbits. The actual semiclassical 
amplitudes satisfy ( |66|) close to the bifurcation, but the agreement is satisfactory in 
the immediate neighbourhood of the bifurcation only. 

For p{r'^) in equation ( |5^ ) we choose a first-order polynomial 

p{r^) ^ Po + Pi{r^ - a) (67) 

for the amplitude function. The coefficients po and pi are chosen to satisfy (|6^: 

P° = AT^T-," , ■^ ' Pi " (aA^^ + Aion ) • (68) 

V27r7r(l -I- 1) 27ria \ ^2n J 



By (|66|), pi remains finite at the bifurcation. Thus, the uniform approximation ( p6| 
assumes the form 



*(i;) 
with the integrals 



■^""■^ Io + - (aA^^ + A non I -'I 



(1 + i) a \ V2^ 



e'^-- (69) 



^0 = iUT^ j dxdye'^-^-'y^ , = ^ J dx dy {r^ ~ a) e'^^^-^^y^ . (70) 

Jo can be evaluated in terms of the Fresnel integrals | [T^ 

nX nX 

C(x) = / cos {^t^)dt, S{x)= sin {^t^)dt (71) 



to yield 



C = - iS 



(72) 
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Figure 4. Uniform approximation for the generation of the (4,5) non-axial orbit 
from the downhill orhit for the electric field strength F = 10~*. Solid line; 
Uniform approximation (p9|), dashed line: isolated-orbits approximation. 



Ii can be reduced to 

h = 



1 

2i'' 



(73) 



with V = r"^ — a. The integral in ( ^3| ) diverges at t; = cx). It can be regularized by 
adding a small exponential damping factor e^"^^ and letting e — > at the end. This 
procedure yields 

/i = 1 . (74) 

The regularization can be justified by noting that the stationary phase approximation 
to ( [70| ) also yields /i = 1, so that ( |7^ ) has the required asymptotic behaviour. 

The calculation of the uniform approximation ( |56| ) is thus finished. As an 
example, the uniform approximation for the bifurcation of the (4,5) non-axial 
orbit off the downhill orbit is compared to the simple closed-orbit formula in 
figure ^. Obviously, the uniform approximation (69) smooths the divergence of the 
isolated-orbits approximation and, as desired, asymptotically reproduces the simple 
approximation perfectly. 



5. Semiclassical spectra 

A low-resolution semiclassical photo-absorption spectrum can be obtained from the 
closed-orbit sum (^) by including orbits with an orbital period up to a maximum 
Tmax only. In order to resolve individual energy levels, Tmax must be larger than the 




Figure 5. Bifurcations of non-axial orbits off the downhill (+) and uphill (x) 
orbits in a scaled period vs. scaled energy plot. Bifurcation energies and scaled 
periods of the bifurcating orbits are indicated. Dashed lines: single, double and 
triple scaled perturbative Heisenberg time (PSl). 



Heisenberg time Th- A rough estimate for Th can be obtained from perturbation 
theory. To first order in the electric field strength, the energy splitting between two 
adjacent spectral lines with principal quantum number n is p3| AiJp = 3ni^, so that 
the scaled perturbative Heisenberg time is 

Th,p = ^ = ^v/^. (75) 

This estimate is reasonable (although it may not be quantitatively precise) as long 
as different rt-manifolds do not overlap. If they do, the mean level spacing is much 
smaller and the Heisenberg time therefore much larger than given by (|75|). 

Multiples of the scaled perturbative Heisenberg time ( [TSf ) are plotted in figure 
H together with the scaled energies and scaled periods of bifurcating orbits. _For 
low scaled energies, the periods of bifurcating orbits are well approximated by Tjj.p. 
Therefore, there is no parameter range where the closed-orbit sum can be extended 
up to the Heisenberg time without involving bifurcations, so that the uniformization 
of bifurcations must be an essential ingredient to any closed-orbit theory quantization 
of the Stark effect. 

These findings can be confirmed numerically. Figure ^ displays the oscillatory 
part of the photo- absorption cross section as calculated from the closed-orbit sum. All 
spectra presented in this section are non-scaled spectra calculated for the hydrogen 
atom in an electric field F = 10~*=51.4V/cm, corresponding to w = 100 and 
E — 10"'£', with the initial state |ls) and light linearly polarized along the electric 
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Figure 6. Closed-orbit sum without uniform approximations as a function of 
energy at the electric field strength oi F = 10~* with scaled cut-off time (a) 
Tmax = 5, (b) Tmax = 6.5. At E > -5 X 10"*, the latter is dominated by 
bifurcation-induced divergences. 
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field axis. Small discontinuities are introduced in the low-resolution spectra because 
a closed orbit abruptly disappears from the truncated sum when its period increases 
beyond the chosen cut-off time T^ax- The discontinuities can be avoided by choosing 
a smooth cut-off function to gradually switch off the contribution of an orbit when its 
period approaches the cut-off time. The changes in the spectra are hardly visible, and 
as excellent results can be obtained with the simple hard cut-off, no smoothing will 
be used henceforth. ^ 

For the spectrum in figure ^(a), a scaled cut-off time of Tmax = 5, well below the 
perturbative Heisenberg time, was chosen. Consequently, one can distinguish groups 
of levels characterized by a fixed principal quantum number n, but no individual 
spectral lines can be resolved in the plot. The small oscillations are an artefact of the 
closed-orbit sum and cannot be identified with spectral lines. This is clear from the 
observations that the frequency of the oscillation increases with increasing cut-off time 
and that virtually the same oscillations are visible within an n-manifold and between 
manifolds. 

At low energies, different n-manifolds do not overlap, so that spectral regions 
with high oscillator strength density alternate with regions where no spectral lines 
are present. (As the smooth part of the semiclassical spectrum has been omitted, the 
bottom line of the semiclassical spectrum is shifted to negative values of the oscillator 
strength density.) At E ^ —4.7 x 10^^, neighbouring n-manifolds start to overlap. 
In the overlap regions the spectral density is considerably higher than in regions of 
isolated n-manifolds. 

To improve the resolution, the cut-off time must be increased. Figure ^b) displays 
the results for a scaled cut-off time of Tmax = 6.5. At low energies, this is still below 
the perturbative Heisenberg time, and no significant improvement of the resolution 
can be found. At high energies, bifurcations start to occur, and a dense sequence of 
bifurcation-induced divergences covers the semiclassical signal. As the cut-off time 
is further increased, bifurcations occur at ever lower energies and destroy ever larger 
parts of the semiclassical spectrum. It is thus obvious that in its simple form the 
closed-orbit sum is useless for a complete quantization. 

If, on the contrary, the uniform approximations are included, the closed-orbit 
sum can be extended to longer orbits. Figure |^ shows the uniformized closed-orbit 
sum for a scaled cut-off time Tmax = 15. As will also be done in all subsequent 
semiclassical spectra, the uniform approximation was applied for bifurcating orbits 
whose action difference is less than 2tt. If for a given orbit several other orbits 
satisfy this requirement, a more complicated uniform approximation describing several 
bifurcations should be used (see the discussion at the end of section ||). For the time 
being, we resolve this conflict by calculating the simple uniform approximation for the 
two orbits with the smallest action difference and treating all other orbits as isolated. 

Although individual spectral lines can be discerned in figure ^ it would be hard 
to obtain precise values of the energy levels and, in particular, of the associated dipole 
matrix elements from this figure. A more reliable method of extracting the spectral 
information from the semiclassical data is clearly desirable. To this end, we will use 
the method of semiclassical quantization by harmonic inversion 
recently generalized to be applicable in connection with uniform approximations |12| . 
This method requires the cut-off time of the closed orbit sum to be chosen larger than 
twice the Heisenberg time for this purpose. The perturbative Heisenberg time thus 
suggests that a scaled cut-off time of T^ax = 15 will be sufficient as long as different 
n-manifolds do not overlap. Figure presents the energy levels and transition matrix 



,18, la, which we 
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elements obtained from it by the harmonic inversion procedure and compares them to 
the exact quantum results. 

In the region of isolated n-manifolds the quantum and semiclassical spectra can 
indeed be compared line by line to reveal excellent agreement in both the position 
and the intensity of spectral lines. Note, in particular, that the low-resolution 
semiclassical spectrum apparently ascribes considerable amplitudes to spectral lines 
in the middle of an n-manifold, whereas the quantum amplitudes nearly vanish. The 
high-resolution semiclassical spectrum correctly identifies these amplitudes as being 
very small. Furthermore, it can be seen from the magnified low-resolution spectrum of 
the manifold rt = 31 in figure ^a) that the spectral lines, in particular the weak lines 
close to the centre of the manifold, appear asymmetric in shape in the semiclassical 
spectrum: They extend farther to the right than to the left. Nevertheless, the high- 
resolution analysis identifies the lines correctly. 

Figure ^(b) shows the low-resolution spectrum for the manifold n = 31 with a 
signal length of Tmax = 40. It can be seen that the spectral lines close to the centre 
of the manifold become smaller in comparison to the outer lines, thus approaching 
the true semiclassical spectrum, but that the asymmetry of the lines is exacerbated: 
They are clearly saw-tooth shaped, rising steeply on the left and falling off gently to 
the right. 

At E K —4.7 X 10"'', neighbouring n-manifolds overlap for the first time, thus 
doubling the density of spectral lines. It can be seen in figure || that at this energy 
the harmonic inversion of the given semiclassical signal abruptly breaks down. It 
recovers at slightly higher energies, where again only levels of a single n-manifold are 
present. At « —4.5 x 10~^ and E « —4.2 x 10~^, pairs of levels belonging to 
different n-manifolds are so close to being degenerate that they cannot be resolved by 
the harmonic inversion. In the semiclassical spectrum they appear as single lines with 
amplitudes equal to the sum of the two quantum amplitudes. 

Similar effects can also be observed at higher energies. However, as the energy and 
the density of spectral lines are further increased, the harmonic inversion gradually 
ceases to yield meaningful results. In the high-energy region at i? « —3.2 x 10~^, 
a few lines can, somewhat arbitrarily, be identified, whereas most of the lines from 
the quantum spectrum are absent. In this region the cut-off time of the semiclassical 
signal is evidently, and expectedly, too small. 

The harmonic inversion always yields some spurious spectral lines together with 
an error parameter e that can be used to distinguish between true and spurious lines Q . 
In figure ^, only lines with an error parameter e < 6 x 10~^ have been included. This 
threshold value will also be used for all other semiclassical spectra presented in this 
section. The selection of "good" semiclassical eigenvalues is therefore solely based 
on criteria inherent in the semiclassical quantization procedure, no lines are selected 
according to how well they fit the quantum results. As can be seen from figure H, 
some spurious lines pass the selection. They are all characterized by having small 
amplitudes, so that they exert little influence on the semiclassical signal. Only rarely 
does it occur that a true spectral line of considerable amplitude is removed. If the 
region of overlapping n-manifolds, where the signal is too short to resolve the spectral 
lines, is ignored, the only instance of this error present in the above spectrum can be 
found at Ek -5.1 x 10"'*. 

To improve the resolution, the cut-off time needs to be increased. As in the 
spectral region around E « —3.2 x 10~* three (or even four) different n-manifolds 
overlap, the true Heisenberg time can be expected to be close to three times its 
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Figure 8. High-resolution semiclassical (sc) and quantum (qm, inverted) photo- 
absorption spectrum. The scaled cut-off time for the semiclassical spectrum is 

^max — 15. 
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perturbative value. Figure || suggests choosing a scaled cut-off time of T^ax = 40. 
Results obtained with this semiclassical signal are shown in figure In the low-energy 
spectral range, the amplitudes of the spurious lines have diminished. Thus, even in 
this range the longer signal yields a better semiclassical spectrum. In the region of 
overlapping n-manifolds around E « —3.8 x lO""*, all spectral lines are well resolved. 
At -E « —3.2 X 10~^, groups of three closely spaced lines can be identified in the 
quantum spectrum. In many cases, all three of them are resolved in the semiclassical 
spectrum. In a few cases, nearly degenerate lines appear in the semiclassical spectrum 
as a single line whose amplitude is the sum of the two quantum amplitudes. 

The most difficult spectral region appears at E ^ —4.2 x 10^^. In this region, 
spectral lines belonging to neighbouring n-manifold are so close to being degenerate 
that they are hard to distinguish even in the quantum spectrum. This degeneracy is 
somewhat accidental, as it will be lifted as the electric field strength is varied, but 
it nevertheless poses a particular challenge to the harmonic inversion. Even with the 
signal length of Tmax = 40, the degenerate levels cannot be resolved semiclassically. 
Their resolution would presumably require a significantly longer semiclassical signal. 
As an alternative, the harmonic inversion of cross-correlated signals has proven 
powerful in resolving nearly degenerate levels |2^, It can be combined with the 
novel quantization procedure for non-scaling systems in an obvious way and can be 
expected to considerably reduce the signal length required to identify the unresolved 
spectral lines. 

Results obtained with Tmax = 30 for lower energies are shown in figure ^ In this 
region all n-manifolds are well isolated. As the splitting of levels within a manifold 
is smaller than at higher energies, the perturbative Heisenberg time is Th,p — 9.4 at 
E — — lO^"'. Therefore, the chosen signal length should easily suffice to resolve the 
spectral lines. This is indeed achieved, the quantum and semiclassical spectra can 
be compared line by line. (Most of the small spurious lines between the manifolds 
can be removed if a stricter quality criterion for the semiclassical lines is applied. In 
figure |ll|, the same threshold e < 6 x 10 * as in previous figures is used.) However, the 
agreement between the spectra is not as perfect as was found at higher energies. In 
particular, considerable contributions to the semiclassical spectrum arise at the centres 
of the n-manifolds at low energies. This decrease in quality should be expected when 
the semiclassical approximation is applied to low-lying states. This is a fundamental 
limitation of semiclassical methods. In this case, however, a more subtle difficulty 
should be considered: Any non-axial orbits suffers two bifurcations, viz. its generations 
from the downhill orbit and its destruction at the uphill orbit. At low energies, 
the distance between the two bifurcations gets small, so that eventually they can no 
longer be regarded as isolated. At this point, the uniform approximation constructed 
in section ^ must fail. It must then be replaced with a uniform approximation 
describing both bifurcations collectively. Such an approximation can be constructed, 
but, contrary to all previous uniform approximations, it must use a spherical rather 
than Cartesian configuration space. This will be discussed in detail in a forthcoming 
pubhcation [Q. 

An even harder uniformization problem arises close to the Stark saddle energy. 
As this energy is approached from below, the downhill orbits undergo an infinite 
sequence of bifurcations in a finite energy interval. Similar bifurcation cascades have 
been observed in the hydrogen atom in a magnetic field and in Henon-Heiles type 
systems . Available techniques for the construction of uniform approximations can 
only take a finite number of bifurcating orbits into account, so that they do not allow 



The hydrogen atom in an electric field: Closed-orbit theory with bifurcating orbits 25 



CM 

c 
Q 



CO 
o 



15 
10 
5 

-5 
-10 
■15 

10 



n 1 1 1 1 1 1 r 



(sc) 



(qm) 



_l I I I I I L. 



-9.6 -9.4 -9.2 -9 -8.8 -8.6 -8.4 -8.2 -8 -7.8 

lO^E 



c 
Q 

CO 



CO 

o 







-5 



-10 



(sc) 



(qm) 



-7.6 -7.4 -7.2 -7 -6.8 -6.6 -6.4 -6.2 -6 -5.8 



lO^E 



Figure 11. High-resolution semiclassical (sc) and quantum (qm, inverted) photo- 
absorption spectrum. The scaled cut-off time for the semiclassical spectrum is 

^max — 30. 
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the uniformization of an infinite cascade. For this reason, the spectral region around 
the Stark saddle energy was excluded from the calculations of the present paper. 

6. Conclusion 

In this paper we have solved the long standing problem of semiclassical quantization of 
the hydrogen atom in an electric field based on closed-orbit theory. We have presented 
a detailed study of both the classical dynamics of the hydrogen atom in an electric field 
and its semiclassical description in the framework of closed-orbit theory. The classical 
closed orbits were described, explicit formulae for the important orbital parameters 
were given, and the pattern of closed-orbit bifurcations was derived analytically. We 
constructed a particularly convenient form of the uniform semiclassical approximation 
for the closed-orbit bifurcations in the Stark system. Contrary to earlier versions, our 
uniform approximation can be computed from parameters characterizing the isolated 
orbits only, so that it is as easy to apply as the simple isolated-orbits approximation 
and can immediately be included into a closed-orbit summation. 

We demonstrated that a semiclassical quantization of the hydrogen atom in 
an electric field by means of closed-orbit theory is impossible unless uniform 
approximations are included. By means of a generalized method of semiclassical 
quantization by harmonic inversion introduced recently , we were able to calculate 
high-quality semiclassical photo-absorption spectra from the uniformized closed-orbit 
sum. We thus proved that our novel method is indeed well suited to the extraction of 
high-resolution semiclassical spectra from low-resolution spectra of arbitrary systems, 
whether they are classically integrable, chaotic, or even mixed regular-chaotic. It thus 
paves the way to a semiclassical quantization of the hydrogen atom in a magnetic field 
or in crossed electric and magnetic fields. Further work on these more complicated, 
classically chaotic systems is in progress. 

The present uniform approximation is applicable only as long as individual 
bifurcations can be regarded as isolated. It will fail if they are close. For this reason, 
the description of the spectral region close to the Stark saddle energy remains an open 
problem. At this energy, the downhill orbits undergo an infinite cascade of bifurcations 
in a finite energy interval. The uniformization of this kind of cascade is far beyond 
present-day techniques for the construction of uniform approximations. It presents a 
challenge worth-while for future research, both as a problem in its own right and in 
view of its potential applications: The bifurcation cascade occurs at the continuum 
threshold, where the transition from bound states to a continuum with embedded 
resonances takes place |2). The semiclassical description of this transition region is 
a task of special importance. It will require the uniformization of the bifurcation 
cascades. 
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